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We give an elementary discussion of parton saturation and its descrip- 
tion by the effective theory of the Color Glass Condensate. We report on 
progress in calculating multi-gluon correlators. The latter are necessary for 
(— | many phenomenological applications, upon some of which we briefly touch. 

Q_i 1. From partons to the Color Glass Condensate 

<D 

A hadron is a complex object in its rest frame as valence quarks are 
accompanied by a sea of hadronic and vacuum fluctuations. Both types of 
y—i such virtual fluctuations are of non-perturbative nature and have the same 

lifetime AtRp ~ 1/Aqcd, since Aqcd is the only available energy scale. The 
picture changes in the infinite momentum frame; the lifetime of the vacuum 
ones remains the same (the vacuum is Lorentz invariant), but the hadronic 
C*~) ones are time-dilated, that is AtiMF ~ 7/ Aqcd with 7 > 1, and the two 

types are disentangled [l]. Probing the hadronic fluctuations is most easily 
done by DIS of a lepton off the hadron. The collision time is estimated as 
Ai C oll ~ 2xP/Q 2 . Thus, quarks with transverse momenta k± <C Q, will be 
"seen" by the virtual photon, since Aifl uct ~ 2xP/k\ > Ai co n. 

The density of partons depends on their momenta and is determined by 
the QCD dynamics. The probability for a parton in representation R with 
!h momenta (0±,p z ) to emit a soft gluon with (k±, k z = xp z ) at small-x is 

dP = c fl a # )d ^ dx . (1) 

7T J_ *^ 

The emission of soft and collinear gluons is favored, since it goes with a large 
logarithm which may even overcome the smallness of the QCD coupling. 
Such radiative processes need to be resummed to all orders and lead to the 
DGLAP and BFKL evolution equations, in InQ 2 and lnl/x respectively, 
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and the one to be used depends on the kinematics. The last factor in Eq. ([T]) 
is present because the emitted parton is a gluon, hence the small- x tail of 
the wavefunction is dominated by gluons. For a BFKL cascade generated 
by a single valence quark one finds the unintegrated gluon distribution, 
i.e. the number of gluons per unit rapidity Y = lnl/x and for a given k±, 
to increase exponentially with Y. 

DGLAP and BFKL evolutions are linear and emissions from the source 
partons take place independently but they lead to totally different outcomes, 
cf. Fig. [lja. Increasing Q 2 , DGLAP produces more partons (for x < 0.1), 
but they occupy a transverse area l/Q 2 and thus the partonic system be- 
comes more dilute. It is self-consistent in the sense that linear dynamics 
remains valid in the course of evolution. On the contrary, increasing In 1/x, 
BFKL produces partons typically of the same size which eventually overlap. 
Beyond that point the source partons emit coherently, the BFKL equation 
becomes inadequate and has to be supplemented by non-linear terms which 
lead to parton saturation. A simple criterion to define saturation is via the 
gluon occupation number, the number of gluons at a given x multiplied by 
the area each gluon fills up and divided by the transverse hadron size. We 
define the saturation momentum as the line along which this occupation fac- 
tor is constant and of order l/a s , namely xg(x, Q 2 (x))/[Q 2 (x)R 2 ] ~ l/ce s , 
with R the hadron size and xg the integrated gluon distribution. Roughly, 
Q 2 is increasing as a small power of 1/x {2|j4j, reflecting the need to stop 
the exponential growth in Y discussed earlier and for a large nucleus ispro- 
portional to A ' , since xg grows with A and R A with A 1 / 3 . Cf. Fig. njb. 




A modern effective theory encompassing parton saturation is the Color 
Glass Condensate (CGC) [l]. Color since we are dealing with QCD, glass 
since there is a separation in times scales-source partons live longer than 
emitted ones-, and condensate since occupation numbers can be as high as 
l/a s . It describes the hadron in terms of a classical field A which reaches 
values of order 1/g. Parton distributions are replaced with a classical proba- 
bility distribution VFy[*4.] encoding higher-point correlations. The evolution 
of this functional is quantum and obtained by summing all orders in aY, 
with a = a s N c /ir, like in BFKL, and all orders in the field A. One has (5}|7| 



with 0[A] a generic observable and H the JIMWLK Hamiltonian. 



The propagation of a projectile parton through the CGC is eikonal. With 
the parton at transverse position x± and moving along x~ , the interaction 





2. Particle production in pA collisons 



The Color Glass Condensate and some applications 



3 




with the field Aa = b~^ + a a in covariant gauge is given by the Wilson line 
V x = Pexp hg f dx~a a (x~, x±)t a ] . We want to study particle production 
in hA collisions, with h a dilute projectile hadron and A a dense target, like 
a large nucleus. The partonic process is qA — > qX, cf. Fig. [2ja, and, in the 
proton fragmentation region, a large- x quark from the proton interacts with 
the soft components of the nucleus via the Wilson line. The cross section is 



dy d 2 p A 



xif q (x hP 2 ± ) J d 2 r ± e- ir ^ (S X1X2 ) Y . (3) 



Here p± and y are the transverse momentum and rapidity of the produced 
quark, x\f q is the quark distribution in the proton, r± = x\ — X2, S XlX2 = 
(1/N C ) tr(y xi V X2 ) with the Wilson lines in the fundamental representation 
is the "dipole" S'-matrix and the average is to be taken with Wy [a] of the 
nucleus. The fractions x\ = p±e y /\/s and X2 = e~ Y = p±e~ y / \fs <C 1 are 
determined from the kinematics in the Born approximation. The ratio 

_L_ dN h /dr]d 2 p ± \ pA 
UpA A 1 / 3 dN h /dr)d* P± \ pp [V 

should be suppressed compared to unity since a big nucleus is more saturated 
than a proton. RpA is measurable and fits to RHIC data from dAu collisions 
are successful and predictions for the pPb ones at LHC are available [8j. 

Similarly we can study di-hadron correlations in hA collisions in the 
forward region. One expects a suppression of the azimuthal correlation of 
the two hadrons when their transverse momenta are of order 0(Q S ), and 
in fact this has already been observed and explained by CGC at RHIC 
J9j[l0]. Studying the process at the partonic level one can derive the anal- 



ogous to Eq. g cross section for the production of, e.g., a qg pair. It is 
straightforward to see, from squaring Fig. [2]b, that it is necessary to in- 
troduce higher-point correlators involving the product of the quadrupolc 
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(a) (b) (c) 

Fig. 2. (a) Quark production and (b), (c) quark-gluon production in pA collisions. 

Q Xl x 2 x 3 x4 = (1/N C ) tr(V r i 1 T4 2 V r i 3 V a;4 ) times a dipole. The average is to be 
done at a large- 1" determined by the kinematics of the process. 



3. Multi-gluon correlators in JIMWLK 

The JIMWLK Hamiltonian determining the QCD evolution of Wy[a] 
and hence of the (gauge invariant) correlators is 

with M. uvz = (u— v) 2 /[{u— z) 2 (z— v) 2 ] the dipole kernel, the tilde standing 
for the adjoint representation and the functional derivatives acting on the 
upper and lower end-points of the Wilson lines and V respectively. 

One method to calculate the correlators is to reformulate Eq. ^ into a 
Langevin equation [IT] for a Wilson line and solve on a lattice. The other 
is to construct the evolution equations for the correlators (0)y of interest 
using Eqs. ^ and ([5]). The structure of the real terms (the last two) of the 
Hamiltonian is such that it leads to a hierarchy of equations [12] ; H reSL \ acting 
on n Wilson lines, can lead to n+2 of them. This second method is sufficient, 
if, for instance, one wants to calculate single inclusive gluon production in 
pA collisions, cf. Eq. ([3]). JIMWLK leads to the dipole equation 



0(S XiX2 )y _ a_ 
dY ~ 2tt 



/ J^x\X2Z {Sxiz£*zx2 Sxxx^Y j (6) 
J z 



which is interpreted in terms of projectile evolution as shown in Figj3la at 
large- N c . One factorizes (SS)y = (S)y(S)y to get the BK equation 12ll3 , 
a closed equation whose solution is known semi-analytically and numerically. 

Less inclusive quantities require the knowledge of higher-point correla- 
tors. The quadrupole equation [l4] is more involved and two representative 
diagrams are shown in Figs.[3}b and[3jc. Even at large- N c , where it becomes 
a closed inhomogeneous and linear equation (requiring the input of (<S)y), 
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Fig. 3. (a) Real and virtual emission of a qq pair from a dipole. (b) and (c) Similarly 
from a quadrupole. The dashed line denotes the interaction with the target. 



the large number of transverse variables and the non-locality in transverse 
space seem to be prohibitive for the possibility of a numerical solution. 



The implementation of the Langevin approach 15 , lead to the finding 
that the numerical data for the quadrupole are very well-described by a 
Gaussian mean field approximation (MFA); an extrapolation to arbitrary 
Y of the McLerran-Venugopalan (MV) model, the typical initial condition. 
Such a Gaussian Wy[a] involves a single kernel, hence all high-point correla- 
tors are expressed in terms of the 2-point one, although there was no a priori 
reason for this to happen for the, highly non-linear, JIMWLK evolution. 

Still, one can prove that a Gaussian approximation is a quasi-exact so- 
lution to the JIMWLK equation 16 17]. At saturation, real emissions are 



suppressed and the virtual part H vir t (the first two terms) of the Hamilto- 
nian dominates. This is Gaussian, including the second term; those Wilson 
lines transform the "left" derivatives to "right" ones which act on the lower 
and upper end-points of the Wilson lines and V respectively. We inte- 
grate the kernel over z for 1 /Q s <C \u — z\, \v — z\ <C \u — v\, with the lower 
limit imposed by our approximation, while the upper gives the dominant 
logarithmic contribution 21n[(it — v) 2 Ql]. This would lead to a result valid 
only at saturation, but doing the same approximation in Eq. ^ we find this 
logarithm to be related to the logarithmic derivative of the dipole. Thus 

tt __ 1 f dln(S uv ) Y f_5 6_ _6 5_\ 

G 4g*C F J uv dY { 5al u 5a" Lv + 5c& u 5a a Rv )' U 

This is a Gaussian Hamiltonian and the kernel can be most easily deter- 
mined from the BK equation. Eq. is correct, at finite- N c , at saturation 
by construction and in the dilute limit as can be inspected. 

Hq generates evolution equations for multi-gluon correlators which are 
local in the transverse plane, that is, they are ordinary differential equations 
in Y with Y-dependent coefficients. A further "separability" property of 
the kernel in Eq. Q leads to correlators which are local function in Y . For 
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Fig. 4. (a) (5 6 )y = (^V c 2 /(^ c 2 " IXQ1234S41 " (1/N^)S 32 ) Y with ^3 = x x and 
a; 4 = x 2 . Red: JIMWLK from F = to 5.18 15]. Blue: MFA analytic result for 
N c = 3 [17]. At y = JIMWLK and MFA coincide because of the initial condition. 
True at any Y in the dilute and in the dense region. A tiny difference, stabilizing 



soon, occurs in the transition region, (b) Scaling of the gluon spectrum 18 . 



instance, the quadrupole at large- iV c with MV model initial condition reads 

In [{Si 2 )y{S34)y/{Si3)y{S24}y] 



(Ql234)Y 



+ 



In [{S 12 )y{S34:)y/{Su)y{S23)y] 
In [{Su)y(S23)y/ (Si 3 )y {S24)y\ 
In [{Su)y{S23)y / (Si2)y(S34)y] 



(Si2)y(S34)y 
(Su)y(S23)y 



(8) 



(with i = Xi). It obeys a mirror symmetry (Qi234)y = (Qi432)y> which 
holds at finite-A^ c and beyond the MFA and is a result of time-reversal 
symmetry, where time stands for x~ 17 . It is preserved by JIMWLK due to 



the presence of both left and right derivatives and suggests that the hadron 
expands symmetrically in the x~ direction. At the level of the MFA only, 
it is also symmetric under the charge conjugation (Qi23i)Y = (Q234i)y- 

At finite-Af c operators mix, e.g. the phenomenologically interesting 6- 
point operator QS mixes with two more and the equation from the matrix 
diagonalization is cubic. As shown in Fig. [4} a for a simple configuration 
the agreement of the analytic and numerical results is excellent. Analytic 
expressions for multi-gluon correlators are necessary to reduce the numerical 
cost for Fourier transforming expressions like Eq. g (and more complicated 
ones, like in the double inclusive case) and obtain the desired cross sections. 

4. More properties and applications 

One expects quantities to scale with Q s [l9j[20]. This property, first 
observed in DIS I21J, extends to other quantities, like gluon production at 
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early times in AA collisions 18 , cf. Fig. [4jb. In such collisions, due to final 
state interactions, the CGC cannot describe spectra of produced hadrons, 
but it can give more inclusive quantities, e.g. total multiplicities |1|. 
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